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1.  INTRODUCTION 


According  to  the  contract,  the  goals  of  this  work  are: 

•  search  in  the  literature  for  and  selection  of  sources  of  chaotic  oscillations  in  radio  fre¬ 
quency  band  with  known  models  and  which  potentially  could  be  used  for  ranging  and  de¬ 
tection  problems. 

•  study  of  mathematical  models  of  chaotic  sources,  identification  of  modes  of  generation  of 
low-dimensional  chaos  and  study  of  modes  of  chaotic  generation. 

•  analysis  of  a  possibility  of  symbolic  description  for  systems  under  consideration. 

Symbolic  description  means  a  representation  of  a  continuous  trajectory  of  a  chaotic  generator 
by  means  of  a  finite  alphabet.  At  the  moment,  theoretically  rigorous  construction  of  symbolic 
description  is  possible  only  for  a  limited  number  of  discrete  time  systems. 

There  are  two  methods  for  constructing  a  symbolic  description  known  from  literature: 

•  construction  of  quasi-one-dimensional  maps  from  a  flow  system, 

•  determination  of  a  set  of  unstable  periodic  orbits  (UPO)  of  a  chaotic  attractor. 

Both  methods  imply  the  procedure  of  embedding  of  an  observed  trajectory  into  a  multidimen¬ 
sional  (preferably,  three-dimensional)  phase  space.  Then  a  one-dimensional-one-dimensional 
map  is  constructed  or  a  topology  of  the  attractor  by  means  of  UPO  is  determined.  With  a  suc¬ 
cessful  embedding  and  determining  of  a  UPO  set,  the  known  structure  of  the  attractor  facili¬ 
tates  determining  the  size  of  an  alphabet  necessary  for  symbolic  encoding  of  the  trajectory. 
Finally,  UPO  set  enables  to  determine  a  generating  partition. 

This  is  why  systems  with  low-dimensional  chaos  are  considered  in  the  literature  in  detail.  For 
such  systems  it  is  possible  to  construct  quasi-one-dimensional  succession  sequence,  which 
detennines  generating  paitition  and  enables  to  solve  some  other  problems  concerning  the 
structure  of  attractor.  In  this  paper  attention  is  paid  to  construction  of  quasi-one-dimensional 
maps,  which  might  define  a  generating  partition. 

Such  a  sequence  could  facilitate  comparison  of  signals,  decision  about  identity  of  signals  (fel¬ 
low/stranger),  synchronization  of  two  chaotic  sources.  In  some  cases  symbolic  description 
could  facilitate  synchronization  in  the  presence  of  additive  noise  or  distortion  of  signal  in  a 
channel. 

In  section  2  a  review  is  given  of  five  models  of  microwave  chaotic  generators:  klystron  gen¬ 
erator  with  delay,  generator  based  on  traveling-wave  tube,  based  on  backward-wave  tube. 
Pierce  diode  and  microtriode.  Two  models  have  been  selected  for  further  investigation:  gen¬ 
erator  based  on  traveling- wave  tube  and  generator  based  on  microtriode. 

Section  3  describes  features  of  such  models,  their  characteristics  and  numerical  methods  used 
for  studying  these  models. 

Existence  of  low-dimensional  chaotic  modes  in  selected  models  of  generators  is  shown  in  sec¬ 
tion  4.  '  ■  - 

The  description  of  trajectories  of  models  of  generators,  selected  in  section  2,  by  means  of  bi¬ 
nary  symbolic  sequences,  is  investigated  in  section  5. 


^  the  symbolic  description  is  com  Jcted  and  properties  of 

obtained  symbolic  sequences  are  studied.  Methods  for  constructing  symbolic  binary  se¬ 
quences  closely  connected  with  imtial  chaotic  trajectory  are  indicated  and  properties  of  such 

sequences  are  considered  for  the  system  based  on  traveling-wave  tube.  ^ 


2.  REVIEW  OF  GENERATORS 

tai'  of  generators  of  chaotic  oscillations  in  microwave 

g^erators  wWcb^b  relatively  simple  mathematical  models  of 

generators  which  have  been  shown  to  have  chaotic  modes  are  included  in  this  review.  This 

of  nei&er  defied  dpcription  of  physical  organization  of  devices  nor  principles 

tion^  detailed  matenal  can  be  found  in  literature,  cited  in  the  body  of  this  sec- 

or  physical  experiments  cited  in  the  text  of  this  section  are  taken  from 


2.1  Klystron  generator 

The  scheme  of  klystron  generator  with  delayed 
feedback  is  presented  in  Fig.  1.  Electron  beam 
(shown  witli  arrows)  is  velocity-modulated  in  the 
first  resonator.  Energy  is  taken  off  the  beam  in  the 
second  resonator.  Input  and  output  resonators  are 
connected  with  the  delayed  feedback  line.  The 
input  resonator  is  excited,  by  the  signal  of  output 
resonator  coming  from  feedback  line,  which  is 
characterized  by  modulation  depth  in  and  delay 
time  o  t .  GUQiut  resonator-  is  excited  by  the  cur¬ 
rent  of  the  electron  beam.  If  voltages  at  cavity 
gaps  of  input  and  output  resonators  are  repre¬ 
sented  as  =  Re(y4, 2  (r)exp(ra)Qr)),  the  equa¬ 

tions  for  excitation  of  resonators  become  [1]: 

dAj  cOa  wcDa 

2^' ^ - 50exp(-zcOo60, (la) 
dA.  co„ 

^  +  =  Po-^A(^.Oexp(-zOo//V(,),  (lb) 

where  /  -  is  the  distance  between  the  resonators, 
K  -  wave  resistance  of  output  resonator,  Q,,  Q2 
-  are  loaded  quality  factors  and 


Fig.  1.  Scheme  of  klystron  generator  with 
delay:  1  -  cathode,  2  —  collector,  3  —  elecr 
tron  beam,  4,5  -  input  and  output  resona¬ 
tors,  6  -  accelerating  electrode,  7  -  feed¬ 
back. 


A  0  -  -2//o^,  (IA(f-^/Vf,)IGg/2Fg)  exp(i(p(t  -I/vg)) 


(2) 


Fig  2.  Waveforms,  phase  portrait  and  power  spectra  density  for  two  modes  of  klystron  gen¬ 
erator. 

s  the  first  harmonic  of  the  current  (second  resonator  is  usually  a  part  of  a  high-Q  empty 
esonator,  so  only  the  first  harmonic  of  the  current  exciting  the  second  resonator  is  consid¬ 
ered). 

/,  is  Bessel  function  of  first  order,  (p  =  Arg(A,) ,  0o  =  0l)o7/V(,  -  is  undisturbed  transit  angle 
phase  incursion). 

A^ith  dimensionless  amplitudes  of  oscillations 

F;(T)  =  4(T)eo/2Fo  '  (3a) 

=  (3b) 

vhere  T  =  cogf /\|/ and  vp  =  ooj (5  r  +  /  /  Vo) ,  then 


^  +  y,F;  =7,^2 

dF, 


di 


+  Y2F2=2aJ.(|Fi(T-l)|)-4^ 


(4) 


\F,{x-l)\ 


lere  -  are  damping  parameters,,  a  =mEV|;9oexp[-/(\)/  +  7r/2)],  z  =  KIqI2V^ 

2  amplification  parameters. 

case  of  low-Q  of  input  resonator  (y,  » 1 )  tlie  system  is  reduced  to  a  first-order  equation; 


5 


~ + y^F  =  2aJ,  u'F(t  - 1)  |) 


(5) 


th^?  w  following  values  of  the  system  parameters 

resoltors  f-in^  r  ''  with  velocity  of  light,  distanL  between 

nators  /  - 10  cm,  fundamental  frequency  of  resonators  is  /„  =  3  GHz,  loaded  quality  fac¬ 
tors  5i_2  -300,  wave. resistance  ^  =  100  Ohm,  accelerating  voltage  Fo=2500  V  current 
/o=50,mA.  -  ’ 

Values  of  dimensionless  parameters  are  \|/-i-Tc/2«1807t,  Oq  «207i,'y,2  «1,  s»10'^ 

K^ttaTini  results  of  numerical  and  physical  experiments  (integration  by  Runge- 

Kutta  me^od  of  the  fourth  order)  (see  Fig.2).  Numerical  experiment  demomtrates  fransition 
to  chaos  throng  senes  of  penod  doubling  bifurcations.  The  experiment  corroborates  also  the 
presence  of  a  chaotic  mode.  But  since  the  experiment  was  carried  out  with  five-resonator  klv 
stron,  agreement  with  numerical  experiment  is  only  qualitative. 

The  results  of  numerical  experiment  qre  presented  in  greater  detail  in  [2]  . 


2.2.  Traveling-wave  tube  (TWT) 

Scheme  of  generator  with  delayed  feedback  based  on  traveling-wave  tube  is  depicted  in  Fig.3 


|“  Delay  feedback  — 

^  Output 

—  >ii 

‘h- 

- 1 - - 

—{——4 

3 

- - 

.  Slowing  structure 

4 

Fig.  3.  Scheme  of  TWT-generator:  1  -  cathode,  2  -  anode,  3  -  accelerating 

grid,  4  -  electron  beam 


Advantage  of  TWT  as  compared  to  klystrons  is  amplification  in  a  broad  frequency  band. 
Results  of  experiments  on  searching  for  chaotic  modes  in  a  TWT-generator  are  described. 


First,  generator  was  considered  in  a  general  form  [5,6],  depicted  in  Fig.  4  where  F 
nonlinear  amplifier,  T  -  is  a  delay  and  0)  -  is  a  filter.  ■ 


IS  a 


The  equation  for  oscillations  of  such  generator  is  a  functional  equation: 
^(t)  =  <^F(X(t-T)), 


(6) 


where  <|)j'(?)=  _[^5(x)5-(r-  .,ic,  and  ^(t)  -  is  an 

impulse  response  characteristics  of  the  filter. 

For  narrow-band  filter  and  for  r  =  0  this  equation 
corresponds  to  reduced  differential  equation 

X'  +  zX=^zF{X)  (7) 

or,  in  the  presence  of  delay 

X’  +  zX^zF[X{t-T)],  (8) 

where  s  « 1 . 


Fig.  4.  Block-scheme  of  TWT  gen¬ 
erator 


This^  equation  is  studied  with  cubic  nonlinearity  F  =  yQX(}.-X^)  or  with  nonlinearity  of  the 

type  F{X)  =  aX exp(-(Z  -  p)^ ) ,  which,  according  to  [5],  describes  the  nonlinearity  of  f WT 
in  the  best  way. 

The  results  of  numerical  simulation  of  the  system  are  described  in  [3].  Simulation  was  carried 
out  under  the  assumption  that  the  field  at  the  input  area  of  TWT  is  defined,  phase  wave  veloc¬ 
ity  is  constant  and  is  independent  of  signal  frequency.  For  the  input  area  of  tube  an  approxi¬ 
mation  of  specified  current  is  used. 

Since  the  input  signal  fades  completely 
in  the  local  absorbent  and  the  output 
signal  is  week  and  does  not  influence 
beam,  the  effect  of  volume  charge  is 
neglected. 


based  of 
tube  (carci- 


2.3.  System 
backward^wave 
notron) 

Traveling- wave  tube  is  a  broadband 
amplifier,  but  it  can  operate  as  a 
generator  only  in  a  relatively  narrow 
frequency  band.  As  a  broadband 
g!mjerahaB:;kwal34oka»affdt^^  tufeealljs 
used.  Its  scheme  is  depicted  in  Fig.  .5. 
In  such  generator  amplification  is  pro¬ 
vided  by  interaction  of  backward  wave 
with  electron  beam,  which  is  used  to 
transfer  oscillation  energy  from  the  in¬ 
put  of  decelerating  system  to  the  out¬ 
put,  where  backward  wave  occurs. 
Thus  electron  beam  is  used  simultane¬ 
ously  as  feedback  regulated  by  phase 
(to  change  the  accelerating  voltage), 
which  provides  tuning  of  generation 
frequency  in  a  broad  frequency  range. 

Dynamics  based  on  ,  backward- wave 
tube  is- described  in  [7],  where  there  are 
also  references  to  earlier  works. 


S-'J 

/  ^  J 

tf 

pF 

Fig.'S  Scheme  of  generator  based  on  back¬ 
ward-wave  tube:  1  -  ca^ode,  2  -  accelerat¬ 
ing  electrode,  3  —  periodic  decelerating  sys¬ 
tem,  4  -  collector,  5  -  energy  output,  6  - 
matched  load 


■XO  -20  -to  0.0  to  zo  to 


Fig.  6.  Waveform  and  phase  portrait  of  genera¬ 
tor  based  on  backward- wave  tube. 


■  It  IS  shown  that  generator  b.  j  on  backward-wave  tube  could  be  des  led  by  equation 


dAj  dAj 

— - -  =  aA  A  * 

dx  ’ 


,,  -  9^2 

(l  +  «) - --\ - ^  = 

dz  (9) 

dA,'dA, 
dx  .6^ 

Stem  considerabte  similarity  tvith  finite-dimensional  Lomnz 

2.4.  Pierce  diode 

Pierce  diode  consists  of  two  parallel  infinitely  wide  equipotential  grids  oenetrated  bv 
mon^nergetic  electron  flow.  Charge  density  and  flow  velocity  at  the  input  of  ITStem  ^e 
fixed.  Space  between  grids  are  filled  uniformly  with  motionless  ions.  ^ 

Control  parameter  is  Pierce  parameter  a  =  <o,i/v.  ,  where  .  is  plasma  fte- 

quency.  p.  -  is  charge  density  m  the  flow  at  the  input .  v,  -  flow  velocity  at  the  input.  L  - 
the  distance  between  grids. 

Pierce  has  shown  in  1944  that  uniform  motion  of  the  flow  with  velocity  v(x)  =  v.  becomes 

unstable  at  a  >  ji .  Fluctuations  of  charge  density  in  the  flow  lead  to  redistribution  of  charges 
m  external  cucuit,  connecting  the  grids,  and  to  the  appearance  of  induced  charge  on  lie 
oun  ing  gnds.  As  a  result,  potential  at  the  grids  remains  constant,  and  uniform  field  of  these 

this  perturbation  leais  to  growTof^c^iat:: 
01  cnarge  dei^ity.  Perturbation  evolves  until  a  virtual  cathode  appears  in  the  svstem  -  a  remnn 

bpuSon  fl';“Sr'°“  caflrode  reflects  part  of 

According  to  [8],  Pierce  diode  could  be  used  to  obtain  super-power  microwave  radiation 
Numerical  simuladon  is  based  on  "particle  in  a  cell"  method.  Electron  beam  is  represented  as 
a  body  of  big  particles  (in  one-dimension  case  charged  Sheets),  injected  in  equJ  time  inter 
vals  wtth  consteut  velocity  into  the  interaction  space.  For  every  sheet  an  equation  fs  IXer 

dt^  ^  (10) 

To  calculate  field  strength  and  field  potential  of  volume  charge  and  charge  density  a  uniform 

XLmL  n  of  volume  chargf  in  electostatic  ap- 

proximation  is  defined  by  Poisson  equation:  ^ 

9‘9  2/  ■ 

^  =  a'(p-l) 


O 


Fig.  7.  Dependence  of  dimension  D  on  observation  scale  for  different  dimensions  of  em¬ 
bedding  space  for  two  modes  of  non-regular  oscillations. 


To  calculate  charge  density  a  procedure  of  linear  weighting  of  particles  on  space  grid  is  used. 
Field  strength  is  calculated  by  differentiation  of  obtained  values  of  potential,  The  nurnber  of 

grid  cells  is  128,  and  the  number  of  particles  per  cell  in  unperturbed  state  is  4. 

Analysis  of  oscillatioiis  was  carried  out  by  means  of  observation  of  time  realizations  of  field 
of  volume  charge  E  at  the  output  of  the  system  and  corresponding  power  spectra.  Existence 
of  chaotic  modes  is  shown;  dimension  of  attractor  is  estimated  (see  Fig.7). 
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•  2.5  Oscillator  based  on  vacuum  microtriode 

|tudy  of  chaotic  dynamics  of  vacuum  microtriode  is  de¬ 
scribed  in  [9].  It  is  shown  that  active  microtriode  oscilla¬ 
tor  can  generate  chaotic  oscillations  in  microwave  band. 

Active  oscillator  could  be  constructed  by  introduction  of 
a  feedback  into  the  scheme  of  fnicrotriode  amplifier  in  a 
manner  similar  to  usu^  tube  |enerator.  However,  due  to: 
the  absence  of  saturation  of  nonlinear  chmacteristics, 
self-oscillations  in  such  generator  would  mcfd^e  infi- 
nitely,  leading  to  destruction  of  cathode.  It  is  necessary  to 
introduce  passive  element  with  nonlinear  characteristics 
into  feedback  loop  to.  compensate  the  signal  growth.. 

Scheme  of  proposed  generator  (see  Fig.8)  contains  a  vac¬ 
uum  microtriode  with  nonlinear  oscillation  circuit 
(nonlinearity  is  given  by  resistor),  included  in  a  grid  cir¬ 
cuit,  and  an  inductive  feedback. 

Anode  current  depends  on  voltage  at  the  grid  according  to 
Fauler-Nordheim's  law  of  autoelectronic  emission: 

l^=AF-(aU  +  bfexp[-BO^^^/(F(fU  +  b)],  (12) 

where  4  -  anode  current,  U-  grid  voltage,  A,B  -  almost  constant  parameters,  to  -  work 
function,  F  -  field  amplification,  a  ,b  -  geometric  constants. . 


Fig.  8.  Scheme  of  generator 
based  on  vacuum  microtriode 
with  inductive  feedback  and 
nonlinear  dissipative  element  in 
the  grid  circuit. 


Nonlinear  dissipation  is  giv^  Jy  nonlinear  resistor,  with  characteristic  ) 

/  =  /o[exp(C//£/o)-^l] 

The  generator  equation  in  dimensionless  variables  for  non-autonomm,.  m 


non-autonomous  mode  is  given  by 


^  -K2x  +  2g  +  =  Fsin(j.T) 


is  demonstrated  by  physical  an  experiment  and  numerical  Simula- 

se„.edKfXCri:,^’^  •“  -  P- 


2.6,  Models  selected  for  investigation 

,^ong  five  models  described  above,  capable  of  generating  chaotic  signals  in  microwave 
X  ^  two- ge-’erator  based  on  microtriode  and  generator  based  on  traveling- 

tW  Preferred  models  that  are  studied  to  some  extent  and  at  the  same  time  have  models 
no  require  enormous  computational  resources  for  simulation  of  their  dynamics. 

^  microtriode,  which  is  described  by  non-autonomous 

ordinary  differential  equation  of  second  order.  This  model  is  of  interest  also  because  eenerator 

dimetr  generator  of  chaotic  signal  in  microwave' region  vJth  finite 

dimension^  model  discovered  in  literature.  All  other  above  models  are  descried I  fh^bv 

obtained  which  are  m  good  agreement  with  those  described  in  [9],  The  calculation  of  a  traiec 
tory  could,  be  carr^  out  by  Runge-Kutta  method  of  4-5  order  with  varilrstep  ir^n: 

stramts  are  imposed  on  maximum  size  of  integration  st^.. 

to  °S°^fthe  moM^nf  bf  are  described  by  differential  equations  with  de- 

lay.  i  nese  je  the  model  of  klystron  generator  and  TWT-model.  For  further  investigation  we 
ave  taken  the  system  based  on  TWT,  because  equation  (8)  is  simpler  than  system  (4). 

Lte'e  of  of  TWT-generator  can  be  employed  to 

analyze  of  the  dynamics  of  Wystron  generator  with  delay.  ^  ^ 

Eq^tion  with  delay  (8),  for  fwT-system,  has  infinite  number  degrees  of  fieedom  and  this 
fact  substantidly  comphcates  investigation  of  this  system.  First,  chaotic  modes  of  hivh  di- 
„  might  exist  in  such  system,  and,  as  one  will  see  from  the  results  descril  bdow 
tos  narrows  the  region  of  parameters  useful  to  us.  Second,  calculation  of  equation  (8)  gives 
one  scalar  component,  so  to  study  chaotic  modes  it  is  necessary  to  reconstruct  a  phase  space 

Reconstruction  algorithm  is  described  in  the  next  section.  «  a  pnase  space. 


3.  ALGORITHMS 

for  identification  of  low-dimension. 


By  investigation  of  the  systt  iased  on  microtriode  and  described  b>  ^iiation  (14)  we  do  not 
face  a  problem  of  separating  the  regions  of  high-dimensional  and  low-dimensional  chaos 
Only  simple  chaotic  modes  are  possible  in  the  system.  The  single  aim  of  estimating  of  dimen¬ 
sion  or  the  major  Lyapunov  exponent  is  to  give  evidence  that  the  modes  are  chaotic. 

With  TWT-system  the  situation  is  much  more  complicated.  This  system  has  chaotic  modes  of 
high  dimension  and  isolation  of  “relatively  simple”  modes  of  low  dimension  is  important.  As 
was  already  noted,  calculation  of  equation  (8)  giyes  the  only  scalar  variable,  so,  preliminary 
reconstruction  of  a  phase  space  and  embeddmg  of.  realization  into  this  phase  space  is  neces¬ 
sary  (embedding  algorithm  is  described  below). 

The  first  problem  here  is  that  embedding  space  dimension  for  calculated  signal  is  unknown 
To  estimate  this  dimension  —  the  embedding  space  dimension  —  the  Karhunen-Loeve 
method  could  be  used. 

After  estimation  of  embedding  space  dimension  it  is  possible  to  reconstruct  the  phase  portrait 

and  to  estimate  attractor  dimension  using  correlation  integral  and  to  calculate  Lyapunov  ex¬ 
ponents. 

To  sum  up,  in  this  work  we 

1.  Estimate  embedding  space  dimension  by  the  method  of  Karhunen-Loeve. 

2.  Estimate  attractor  dimension  and  embedding  space  dimension  by  means  of  correlation 
•  integral. 

3.  ■  .Estimate  Lyapunov  dimension  by  Lyapunov  exponents, 

To  illustrate  the  methods  described  and  verify  the  algorithms  employed  we  cite  the  results  for 
Rossler  system 

X  =  -y-z 
^  =  6  -h  (x  —  C)z 

with  parameters  u-  =  0. 1 5 ,  6  =  0.2 ,  c  =  10  corresponding  to  simple  chaotic  mode.  Under 
otlier  values  of  pararrieters  Rossler  system  can  demonstrate  much  more  complicated  dynamics 


3.1.  Estimation  of  embedding  space  dimension  by  Karhunen- 
Loeve  method. 

To  cons^ct  embedding  space  it  is  necessary  to  get  several  coordinates  in  a  multidimensional 
embedding  phase  space  from  one  observed  coordinate  x(t) . 

There  are  several  methods  for  this.  For  instance,  initial  signal  and  its  derivatives  are  consid¬ 
ered  as  coordinates: 

(^(04(0,  ^(0-)’  (16) 

or  time-delayed  values  of  a  single  coordinate  (Takens  method  of  delay  embedding). 

(17) 

Review  of  methods  for  phase  space  reconstruction  by  one  coordinate  can  be  found  in  [12]. 


■  We  use  the  method  of  (  y  embedding.  Its  idea  was  proposed  131  and  mathematical 

CforiM  f “  “  'i'''-  of  justification  is  that  subject  to  some  otSL 

topo  ogical  equivalence  is  shown,  between  the  initial  attractor  of  the  system  and  aftacte  m 
constructed  by  one  observed  coordinate.  :  ^  attractor  re- 

Z  embedding  space  dimension,  let  us  consider  n  values  of  variable 

,1  t  >^to+T-^,o+,(«-i))S  taken  successively  in  time  interval  t,  as  avectorin  n- 

dimensional  phase  space.  Similarly,  X,  = 

Theoretically,  any  value  of  delay  time  x  is  acceptable.  Practically,  too  big  or  too  small 
complicate  analysis  of  reconstructed  attractor.  The  issue  of  optimal 
delay  tme  have  been  discussed  intensively  in  the  literature.  Optimal  delay  time  coSd  be 
taken  as  the  first  zero  of  autocoirelation  function  for  observed  rLizati^ Ts 
mum  of  mutual  information  function,  etc.  Review  of  possible  methods  is  presented  in  [15], 

we  define  delay  time  as  the  first  zero  of  autocorrelation  function  of  the  nroc- 
riod^  Sin  m  ^  rule  of  “3-5  samples  per  pe- 

■  ^  We  also  check  our  results  on  stability  relative  to  sampling  mterval  x . 

Thus,  consider  a  method  for  estimation  of  embedding  space  dimension  described  in  n  71 
and  also  known  in  Je  theory  of  image  processing  as  algorLi 

cording  to  ftis  method,  it  is  neoessaiy  to  cou^uct  u  sAuence  of  vectors  by  xlenl  Stitd 
1,^1,  ■  J  '  where  Nj.-  time  series  length,  n  -  phase  space  dimension, 

to  15"o?20rT??ra  “P^^'ed  embedding  space  dimension 

tsay,  ID  or  2U;.  1  hen  a  matrix  of  trajectory  IS  constructed: 


2=W 


y 

-1/2  ^2 


phas^i^"  The'nuhJ^'^’^fr  '’'■“■'“'i”"®  ot  phase  trajectoiy  on  axes  of  n-dimensional 
phase  space.  Ilie  m^ber  of  linearly  independent  projections,  i.e.,  rank  of  matrix  2  deter- 

m^s  the  number  of  coordinates  required  for  unique  description  of  the  system  state.  .  Ranlc  of 

tnx  2  equals  to  the  number  of  non-zero  eigenvalues  of  matrix  0  =  ZZ  '  . 

Practically,^  due  to  presence  of  noise  or  due  to  the  fact  that  the  system  is  only  approxi- 
mately  fimte-dimensional  eigenvalues  of  matrix  0  will  not  go  to  zero,  but  will  decrease  with 
the  eigenvalue  index  (and  for  white  noise  they  will  decrease  very  slowly).  Thus  it  is  neces- 
sary  to  determine  some  threshold  value  and  to  account  only  for  eigenvalues  ofcov^lc^Z 
tax  0  which  are  bigger  than  the  threshold  value.  u  vanance  ma 

An  example  of  determining  the  dimension  for  Rossler  attractor  is  depicted  in  Fie  9a  Di 
mension  is  estimated  by  the  break  on  the  of  plot  covariance  matrix  eigenvalues.  For  this  ploi 

After  selecting  space  dimension,  it  is  possible  to  choose  new  basis  from  basis  vectors  cor- 
re^ondmg  to  nonzero”  eigenvalues.  Geometrically,  it  means  such  a  rotation  of  axes  of  « - 

ensional  phase  space  when  phase  trajectory  would  have  nonzero  projections  only  on 

it  HXCS. 
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Advantages  of  this  meth  ire  small  requirements  for  computatit  )  resources  and  noise 
resistance.  It  is  shown  in  [19]  that  successive  application  of  this  algorithm  enables  to  get  a 
signal  substantially  cleared  from  noise. 

Besides,  it  is  shown  in  [20]  that  development  of  methods  based  on  decomposition  on  ei¬ 
genvectors  of  covariance  matrix  enables  to  get  estimates  of  attractor  diinension  that  are  bperr 

ating  in  a  broader  region  of  sampling  time  t  than  method  based  on  estimation  of  correlation 
integral,  described  below. 


Fig.  9.  Results  of  estimation  of  embedding  space  dimension  (a)  and  calculation  of  correla¬ 
tion  dimension  (b)  for  Rossler  system  (15). 


3.2.  Estimation  of  embedding  space  dimension  by  correiation 
integral 

The  secoiid  method  for  calculation  of  attractor  dimension  and  estimation  of  embedding  space 
dimension  is  the  method  of  correlation  integrai  pi],  . 

The  idea  of  the  method  is  the  following.  The  number  of  points  in  the  neighborhood  of  radius 
£  grows  with  £  as  8 in  a  space  of  dimension  d . 

To  calculate  dimensions,  in  a  manner  similar  to  that  described  above,  n  -dimensional  vec¬ 
tors  are  taken  from  one-dimensional  time  series  by  Takens  method  and  for  sufficiently  large 
set  of  vectors  (more  than  several  thousands),  a  number  of  vectors  in  the  £ -neighborhood  of 
some  arbitrarily  taken  vector  is  calculated.  An  estimate  d  of  attractor  dimension  is  derived 
from  obtained  dependence  of  the  number  of  vectors  in  £  -neighborhood  on  s . 

Practically,  correlation  dimension  is  calculated  for  different  embedding  space  dimensions 
yi .  Starting  from  certain  n,  calculated  attractor  dimension  does  not  change  (saturates).  This 
value  of  n,  when  saturation  begins,  could  be  considered  as  estimate  of  embedding  space  di¬ 
mension. 

Result  of  dimension  calculation  for  Rossler  attractor  is  given  in  Fig.  9b.  Horizontal  axis 
features  different  space  scales  (different  s).  Vertical  axis  represents  calculated  attractor  di- 


n 


large,  or  ^  j  small  e  error  of  calculation  become.  Jo  large  Estimation  of 

attractor  dimension  IS  taken  at  “average”  8,  where  the  curve  is 

As  can  be  seen  in  Fig.  9b,  Rossler  attractor  dimension  is  slightly  higher  than  two  a«  ran 


3.3.  Calculation  of  Lyapunov  exponents  and  Lyapunov  dimension 

Lyspunov  dimension  is  defined  by 

rfi  =jfc+A±A±::dlii 

l^wl  '  :  ;  ■  ,  (19) 

where  ^  -  Lyapunov  exponents  of  the  system  (in  decreaemg  order)  t  -  maximum  value  sub- 

jectto  >0. 

Supported  by  numerical  evidence,  that  for  “normal”  attractors 
y  p  ov  imension  equals  dimension  calculated  by  correlation  integral. 

moVo  ^  system  based  on  microtriode  one  exponent  is  positive  and  one  is  negative  for  chaotic 
and  thre^''''°^  corresponds  to  attractor  dimension  bSng  between  two 

sufficient  to  calculate  only  the  major  Lyapunov  exponent  It  can  be 

average  ratey^^^^ofthTseiertmL^^^  perturbations  and  calculation  of 

For  the  TWJ-system,  which  is  infinite-dimensional,  such  a  method  is  inapplicable. 

^gonthm  of  calculation  of  major  Lyapunov  exponent  [23]  and  Lyapunov  exponent  spectrum 
[24]  from  observed  time  series  allows  to  overcome  this  difficulty  tLome  extent  ^ 

S'ritbrn'^f  «fjor  Lyapunov  exponent  we  used  our  program.  It  operates  according  to  an  al- 

dresdemmpv  dc/~ttc^cr,/  P  e  ,  cn  is  available  at  htt]3://www.mpipk.s- 

Unfortunately  for  a  TWT-generator  calculation  of  Lyapunov  exponents  requires  considerable 

even  with  ap- 


4.  THE  RESULTS  Jv  ANALYSIS  OF  CHOSEi',  MODELS  OF 

generators 

4.1.  Dynamics  of  the  system  based  on  traveling-wave  tube 


_  ontrol  par^eter  for  a  model  of  TWT-generator  is  parameter  e  featuring  pass-band  of  low- 
trequency  filter.  In  the  beginning,  parameters  a  and  p ,  determining  gain  factor  and  degree  of 

characteristics  were  fixed.  As  can  be  seen  in  bifurcation  diagram  in ' 
undergoes  a  senes  of  period  doubling  bifurcations  as  control  parameter 
decreases  After  the  series  of  bifurcations,  chaotic  modes  follow.  Their  complexity  as  is 

sho^  below,  mcreases  with  decreasing  ofthe  control  parameter.  The  smaller  s  the  more  the 

furcatira^dia^i  richer  the  dynamics  is.  This  can  also  be  seen  in  the  bi- 

shows  results  of  estimation  of  Lyapunov  exponent  spectrum  for  five-dimensional 
been  noted,  calculation  of  Lyapunov  exponent  spectrum  is  a 
omplex  task,  and  results  are  sensitive  to  algonthm’s  details.  So,  the  role  of  Fig  1.0  is  not  to 
pve  exact  values  of  Lyapunov  exponents,  but  rather  to  show  regions  where  one  or  both 
yapimov  exponents  are  positive.  Regions  of  hyperchaos,  where  two  or  more  Lyapunov  ex¬ 
ponents  are  positive,  have  complex  dynamics  and  are  not  considered  here.  We  pay  attention  to 
relatively  simple  low-dimensional  modes,  with  one  positive  Lyapunov  exponent  -  modes  with 

control  parameter  values  between  0.4  and  0.45. 

waveforms,  spectra  and  calculated  embedding  space  dimensions  by 
method  of  Karhunen-Loeve  (nght  in  figures).  Values  of  control  parameter  between  0.75  and 
.  correspond  to  penodic  motion.  As  can  be  seen  fi-om  the  spectrum  of  eigenvalues  of  co- 
vmance  matnx  (result  of  calculation  of  embedding  space  dimension),  two  first  eigenvalues 
de2riSft2amodS°°^  others.  One  may  surmise  that  two-dimensional  space  is  enough  to 

For  more  complex  modes,  for  instance,  for  the  mode  with  8=0.444,  the  difiference  between 
Sfuch  rnST  already  smaller.  In  practice  we  use  three-dimensional  embedding  space 

Figures  16  and  17  show  results  of  calculation  of  correlation  dimension  for  the  above  modes. 
Apparently,  for  s  =  0.1  dimension  exceeds  5  and  saturation  does  not  take  place  even  for  em¬ 
bedding  space  dimension  equal  to  10.  Such  modes  are  not  considered  here. 

For  a  mode  with  s  =  0.4  an  estimate  of  correlation  dimension  is  between  2.5  and  3,  and  satu¬ 
ration  takes  place  when  embedding  space  dimension  equals  three.  In  the  following  section  we 
make  at  attempt  to  construct  symbolic  description  for  relatively  low-dimensional  modes  in 
tbe  range  0.4... 0.45.  As  is  shown  below,  symbolic  description  could  be  constructed  only  in 
rather  narrow  region  of  parameters. 


Fig.  11.  Power  spectrum,  waveform  and  spectrum  of  eigenvalues  of  covariance  matrix 

for  system  (8)  at  a  =  2.05 ,  p  =  1 ,45 ,  s  =  0.75 . 


Fig;  12,  ^  -  Power  spectrum,  waveform  arid  spectrum  of  eigenvalues  of  covariance  matrix 
-for  system  (8)  at  a  =  2.05 ,  p  =  1 .45 ,  s  =  0.5 . 


Fig.  ■  1 0 .  Power  spectrum,  waveform  and  spectrum  of  eigenvalues  of  covariance  matrix 

for  system  (8)  at  a  =  2.05 ,  p  =  1 .45 ,  e  =  0.444 
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Fig.  1 7.  Correlation  dimension  for  system  (8)  at  a  =  2.05 ,  (3  .=  1 .45 ,  s  =  0. 1  (left)  h 

s  =  0.5  (right) 


4.2.  Dynamics  of  the  microtriode  generator. 

Control  parameters  of  generator  based  on  microtriode  (14)  are  V  and  p .  Fig.  1 8  shows  bi¬ 
furcation  diagram  ^d  Fig.  19  shows  the  major  Lyapunov  exponent  as  a  function  of  parame¬ 
ter  F  in  the  range  [0.18. ..0.24]. 


Fig.18.  Bifurcation  diagram  for  system  (14). 


1  n 


AiTipdtude 


Fig.l9.  Diagram  of  major  Lyapunov  exponent  for  system  (14). 

Lyapunov  exponent  indicate  that  a  region  of  parameters 
demonstrates  chaotic  behavior.  In  the  next  chapter  we  construct  svm- 
ynamics  for  this  system  for  parameter  values  V  between  0.222  and  0.233. 

As  an  example,  Fig.  20  shows  the  signal  and  spectrum  of  chaotic  mode  of  system  (14). 


5.  ALGORTIHM  J^OR  CONSTRUCTION  OF  SYMBOLIC 

DESCRIPTION 


In  this  section  we  demonstrate  that  in  some  modes,  continuous  signals  generated  by  discussed 
generators  could  be  associated  with  bin^  symbolic  sequences. 

Let  us  give  basic  definitions  of  the  terms  used,  ■ 

Let  a  map  be  given  by: 

•*^n+i  =  e  G  c  ^20) 

where  R  is  //-dimensional  real  space;  G  -  a  set  in  R^,  where  the  map  (20)  is  defined. 

Any  set  of  m  non-overlapping  regions  P={5,};='",  which  covers  set  G  of  dynamic  system 
states,  is  called  partition 

.  .  (2i) 

Let  us  introduce  a  symbol  /  =  ^Bd,  ieM=  {1.2.3,. ..m]  for  every  region  £,  of  partition  p 
Let  then  denote  a  set  of  all  sequences  of  unbounded  length  X"'=  X^Xj. .  .Xj...,  with  e  M  as 

nZi we  get  a  map  pj3 :  G 

ftp  (ri)  -  X^  ^ \x])  e Bxj,  for  j  >  I,  ^22) 

which  relates  a  sequence  X,"  e  ^  to  every  point  x,  e.  G,  where  Xj  is  a  symbol  generated  at 

time/  Sign  o  means  an  equivalence  between  the  sequence  X"  and  a  sequence  of  sainples 
obtained  after  iteration  of  the  map/xi). 

Accordmg  to  (22)  the  dynamics,  defined  by  a  map/in  the  phase  space,  is  translated  to  a  set  of 
symbolic  sequences  ^  based  on  alphabet  M.  In  such  manner  a  symbolic  dynamic  system  is  a 
left-shift  homomorphism,  d  is  formed:  d(X”),  =  X,,, , (where  X,  -  z-th  symbol  in  X,“ ).  A 
set  of  all  sequences,  generated  by  the  system,  is  called  a  language.  The  study  of  sym¬ 
bolic  dynamic  system  (2^,6),  where  2/,  -  is  a  set  of  all  unbounded  sequences  compatible  with 
I,  is  equivalent  to  the  study  of  the  map,  provided  that  partition  of/is  generating,  i.e.  if  every 
unboundedly  long  symbolic  sequence  corresponds  to  the  unique  initial  condition  xj. 

The  motion  of  dynamic  system  in  continuous  (microscopic)  state  space  is  def.erminp.rl  and 
described  by  equation  (20).  If  ari  initial  state  xi  and  map / are  known,  then  future  evolution  of 
(20)  might  be  calculated.  On  tiie  contrary,  the  motion  of  the  system  on  partition  (macroscopic 
motion)  is  stochastic,  and  trajectories  are  sequences  of  symbols.  If  a  past  trajectory  of  equa¬ 
tion  (20)  on  the  partition  is  known,  it  is  possible  to  predict  its  future  macroscopic  states  only 
in  probabilistic  terms.  Different  states  which  belong  to  the  same  region  Bxj  at  time/  could  get 
into  different  regions  at  timey+l  because  of  stochastic  nature  of  dynamic  system.  State  space 
partition  transforms  a  determined  system  into  a  source  of  information  (messages),  which 

could  be  analyzed  in  terms  of  information  theory.  The  entropy  could  be  derived  for  such 
source  of  messages: 

H;  ~.'gp^x;)\og.P{X’).  (23) 

x!' 
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where  ?(X,)  IS  probabiL ,  of  the  trajectoty  fragment  (word)  X,"  ^1X2... Jsr„.  if  ^  deter¬ 
mines  average  uncertainty  for  prediction  of  a  word  of  length  n  For  the  case  of  known 
symbols,  a  conditional  entropy  of  (..1)  symbols  in  macros^cSjec;:^^:^^^^^^^^  " 

=  if/,  where «=  1,2,. ..  and /zo^  =  .  ^24)’ 

expr^Lr  '’y  "  (20)  f”  tWs  partition  p  is  detennined  by  the 

A' =  lim  Af  -  iim  i//f. 

SSbylS  fofarSie“/S^‘‘“  "  ”  “PP"  ‘™‘ 

/?/:  =  sup /z^ . .  •  . 

P  (26) 

If  h  -  hK,  then  p  IS  a  generating  partition.  An  interesting  property  of  generating  partition  is 

hat  corresponding  map  (22)  is  reversible,  i.e.  if  initial  conditions  do  not  coincide  then  svm- 
bolic  sequences  would  be  different:  from  ^ pp  (r  *)  pp  (x ' *).  . 

In  order  to  map  a  continuous  trajectory  into  a  symbolic  sequence  the  continuous  trajectory  is 
transformed  mto  a  discrete  one  by  introduction  of  Poincare  section  Z  (succession  map)  Then 
c^eoubye  urtersections  of  the  rtajectory  «ith  S  are  considered.  The  seotiorpl^s'ifs™ 
lected  mdividually  for  every  system,  and  not  all  points  of  crossing  of  the  traiectorv  with  Z  are 
used  but  only  those  where  the  trajectoty  crosses  the  plane  in  pattMiS^^ 

Let  illustrate  this  for  the  Rossler  system  (1 5)  with  parameter  values  a  =  0. 1 5 ,  6  =  0.2 ,  c  =  1 0 . 
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Fig.  21.  Phase  portrait  for  Rossler  system  (15)  (left)  and  points  x.Kx..  z,)  of  crossing  of  the 

trajectory  with  the  plane  y  =  0 


^ee-dimensional  phase  porttait  of  the  system  for  such  parameter  values  is  depicted  in  Fig. 
la.  For  this  case,  a  convement  choice  of  the  section  plane  would  be  Z  =  { (;cz)-  y  =  0- 

ine^^/rron“w'“““”‘'  XY  proiection  ani 


In  this  case  crossing  of  the  u  .jectory  with  S  take  place  at ;{:  =  0  and  J  0.  A  map  in  section 
plane  E  introduced  in  this  way  is 

x„+i  =  F(x„)  ^27) 

where  j:  denotes  location  , of  a  vector  on  I,  n  is  discrete  time,  F:  nonlinear  function, 

which  is  unknown  in  general  case.  Map  (27)  for  Rossler  system  is  depicted  in  Fig.  21b' 
where  x„=(x„^„y  Such  dependence  z„(x„)  is  a  result  of  a  good  choice  of  the  section  plane 
and  the  simplicity  of  the  Rossler  system  dynamics  for  given  parameter  values.  For  other  sec¬ 
tion  planes  E  the  dependence  z^ix^)  might  become  more  complex,  and,  in  particular,  lose 

unambiguity  in  the  projection  at  one  or  both  of  the  axes. 

Note  that  thorough  consideration  of  symmetries  and  local  direction  of  the  flow  should  be  done 
separately  for  every  system  under  consideration.  However,  experimentally  observed  low- 
dimensional  attractors  frequently  have  a  simple  (single-loop)  structure  and  choice  of  Poincare 
sections  does  not  present  difficulties. 

In  order  to  get  a  symbolic  encoding  of  the  motion  on  attractors  of  map  F  (27),  a  section  E  is 
di-\dded  by  a  finite  number  6  of  non-overlapping  regions  Bk  with  k  .e  M  =  (0,  1,...  b-l}, 
which  coyer  region,  where  an  asymptotic  behavior  take  place.  Consequently,  a  set  p  =  {Bjc} 
is  a  partition.  Every  trajectory  O  =  {xi,  X2,...x„}  is  associated  with  symbolic  sequence  S' =  5] 
S2  ,  which  consist  of  labels  s,-  of  regions  Bsi ,  which  are  visited  at  time  7=1,2,..  n  In 
turn,  subregions  arc  indexed  by  combinations  of  symbols  ^  si  ...,  and  all  combinations 
begin  from  the  label  Sk  of  parental  set  Bsk-'  thus,  for  instance,  all  points  in  an  element  Bsisi  be¬ 
long  to  Bsi  at  time  /  =  1  and  would  be  mapped  into  regions  Bsi  at  time  i  =  2.  All  points  x  in  Bs 
produce  the  same  symbolic  trajectory  S  under  the  effect  of  ^before  the  expansion  into  differ¬ 
ent  regions  of  p . 

Generating  partition  could  be  constructed  in  a  systematic  way  in  hyperbolic' systems  where 
stable  and  unstable  manifolds  intersect  transversally  in  every  point.  It  is  not  clear,  however 
how  to  identify  a  generating  partition  for  “general”  non-hyperbolic  dynamics,  which  are  char¬ 
acterized  by  tangency  between  stable  and  unstable  manifolds  and  absence  of  exponential  re¬ 
pulsion  in  certain  regions  of  phase  space  (that  corresponds  to  negativity  of  all  Lyapunov  ex¬ 
ponents  calculated  on  finite  part  of  a  trajectory).  Non-hyperbolic  motion  is  typical  for  real- 
world  systems..  - 


Fig.  22.  Quasi-one-dimensional  map  x,.^,  (x, )  for  system  (15)  and  generating  partition 
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to  *e  example  with  Rossi  >stem,  a  generating  partition  is  introd,  ii  through  transition  to 

™s  dependencfis  depicted  in 

As  IS  too™  [28].  a  symbolic  description  can  be  constructed  by  partitioning  of  the  niane  bv 
erh^  hne  crossing  the  maximum  of  the  map,  If  a  viUue  x,  il  fo  die  ^frl  L^Z^ 

me  then  dieyalue  of  the  element  of  the  symbolic  sequence  is  i  (or  «0»).  If  value  x  is 

to  4e  right  from  the  vertical  Une,  then  the  i-th  element  of  the  symbolic  sequence  is  it  '(or 

“responding  to  generating  partition 
I  g  j  p  otted.  Local  minima  above  the  line  correspond  to  symbols  “L”  or  “0”  anri  inr'di 

minima  below  the  line  to  symbols  “i?”  or  “1”.  symoois  >  or  0  and  local 

SSnce  "““I’"™  “  foUow^-  A  discrete  vector 

9  J “tttmuous  chaotic  trajectory  by  means  of  Takensproce- 

ure.  ^  hm,  some  coordinate  x  of  this  sequence  is  considered  and  the  dependence  x,  (x.  1  is 
p  otte  .  this  plot  turns  out  to  be  a  quasi-one-dimensional  unimodal  map  then  it  is  'possibie 

to  define  a  generating  partition:  values  x,  are  associated  with  symbols  i  or  R  («0,.or«l») 

depending  on  whether  the  sample  X,  is  bn  the  left  or  on  the  right  from  the  extremum: 

if  x,<x^ 

As  a  result,  binary  symbolic  sequence  s,  is  associated  with  continuous  chaotic  trajectory 


Fig.  2 j .  Generating  partition  and  X  waveform  of  Rossler  system. 

As  can  be  sho^  if  equations  of  the  system  and  symbolic  sequence  are  known  then  a  chaotic 
jec  oiy  can  be  reconstructed,  in  principle,  with  arbitrarily  high  precision  Algorithm  for  re 
construction  ,s  described  in  detail  in  the  next  section  for  the  sysL  based  onXtSe 


Note  that  in  general  case  holic  sequence  does  not  have  to  be  1  ,iry  for  a  system  v/ith 
rather  complex  dynamics.  But  the  scope  of  this  work  is  deliberately  limited  to  the  simple  case 
of  binary  sequences. 

A  review  pf  methads  for  construction  of  symbolic  dynamics  and  symbolic  partition  can  be 
found  in  [13,  25-27].  In  general  case  construction  of  symbolic  description  is  rather  complex. 
It  can  be  solved  only  in  a  limited  number  of  cases. 

In  section  5.2  we  show  that  for  a  system  based  on  microtriode  a  constructed  sequence  is  sym¬ 
bolic  and  could  be  used  to  reconstruct  (with  arbitrarily  high  precision)  the  initial  continuous 
chaotic  signal.  Thus  we  find  such  a  partition  of  the  phase -space  that  the  obtained  sequence 

(sequence  of  symbols  I,  i?  or  0,1)  corresponds  one-to-one  to  continuous  signal. 

A  similar  result  can  be  obtained  for  generator  based  on  traveling-wave  tube. 


5.1.  Symbolic  description  for  microtriode  system 


Symbohc  description  of  phase  trajectory  for  a  system  based  on  microtriode  is  carried  out 
practically  in  the  same  way  as  for  the  Rossler  system,  but  the  succession  map  is  replaced  by 
the  shift  map  because  the  system  is  not  autonomous. 

A  shift  map  is  constructed  by  means  of  sampling  the  coordinates  over  time  intervals  equal  to 
theperiodof  extemal  force.  Namely,  if  r  is  the  period  then  a  discrete  vector  sequence: 

(Xnyi)  =  (x(t,+iT),y{t,+iT)),  (29) 

is  associated  with  continuous  trajectory  (x(r),y(f))  =  ,  where  is  some  initial  time 

offset.:- 

Fig.  24  shows  a  shift  map  for  the  trajectory  of  microtriode  generator  with  parameters 
p  =  0.l8,  j  =  4.5,  ^  =  10,  go  =0.54,  m  =  0.2,  F  =  0.233. 

As  one  can  see,  die  shift  map  is  a  single-valued  curve,  as  in  the  case  of  Rossler  system. 


Fig.  24b  shows  more  clearly  where  the  points  of  shift  map  lie  on  the  phase  portrait.  Quasi-^ 
One-dimensional  map  x,.^j(x,)  ' (fig.  24d)  is  unimodal  similar  to  the  Rossler  system. 
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Fig.  24.  Shift  map  on  the  phase  portrait  (a);  magnified  shift  map  (b)  and  (c)i  quasi-one- 

airtiensional  map  for  mode  with  F  =  0.233  (d)  in  (14). 

tISeTSSm  Xr  >  quas^ne-dimensional  tmimodal  map,  as 

ror  Kossier  system.  This  sequence  is  associated  with  initio  continuous  chaotic  traiectoiv  A 

predion  this  symbolic  sequence  with  arbitirily  high 


Principal  steps  of  such  reconstruction  are  as  follows, 

symbolic  sequence  ^,,^2,....  of  binary  samples  and  the  system  of  equations  eorrespond- 

mgtoinvestigated  models  of  chaotic  oscillator  be  known.  ^ 

1.  As  IS  shown  in  [29]  for  unimodal  one-dimensional  map  x,„  =/(x,)  sample  values  r 
can  be  reconstructed  from  symbolic  sequence  {s,}  with  arbitrarily  high  precision.  It  can  hi 
reconstructed  owmg  to  knowledge  about  future  values  of  symboUc  sequence  s, .  The  more 
Son  j!  f  “  is  “sed  for  reconstruction,  the  higher  precisi™  of  recon- 

stactionofar,  ts.  Precision  of  reconstruction  depends  eitponentially  on  number  N  of  used 

etements  of  symbolic  sequence.  In  the  case  of  microtriode  the  only  difference  here  is  that  man 

taiLo^*Th”SS“  ““  approximately  based  on  calculation  of  long 

SS'^n  Z  M  “  ‘I  to  calculate  a  map  one  time  with  necessary 

accuracy.  Then  the  map  could  be  reused  repeatedly.  Details  on  this  stage  of  the  akorithm  12 

implementation  and  complications  can  be  found  in  [29]  ■  ^  thm,  its 

2.  Having  obtained  values  x,  it  is  possible  to  obtain  value  y,  by  means  of  Poincare  map 
(shift  map)  (Fig,  24)  which  is  known  numerically. 

3.  Pairs  (x„y)  can  be  used  as  initial  conditions  for  integration  of  the  system  on  time  in¬ 
terval  equal  to  the  period  of  external  force.  Thus,  from  known  (x,.y,)iwe  get  approximate 

trajectory  on  [0,7’].  Known  gives  us  approximate, continuous  trajectory 

on  [r,2r]  and  so  on.  Such  trajectory  can  be  discontinuous  at  moments  of  transition  to  new 

cd^o“ts  ZT  arbitrarily  small  by  more  praise  cal, 


'^/r 


Fig.  25  sho^  the  diffcreni ' '  ietween  the  initial  chaotic  signal  and  1  signal  reconstracted 
ftom  symbolic  sequence  for  the  cases  of  W  =50  and  200  elements  of  symbolic  sequence  s, 
used  for  reconstruction  of  each  point  x^.  ' 


X10' 

(a)  :  (b) 

Fig.  25.  Difference  between  the  initial  chaotic  signal  and  the  signal  reconstructed  from  sym¬ 
bolic  sequence:  (a)  the  length  of  symbolic  sequence  is  iV  =  50  symbols  and  (b)  N  =  200  sym¬ 
bols. 


For  bigger  N  it  is  possible  to  reconstruct  a  tvajQCtoryperfectly'exactly,  i.e.  with  precision  of 
computer  calculations.  • 

Dispersion  a  of  difference  between  the  initial  and  reconstructed  sequences  as  function  of  N 
is.shown  in  Table  1. 


Fig.  26.  Shift  map  and  quasi-one-dimensional  map  for  (14)  with  V  =  0.222 

At  this  point  we  have  shown  an  algorithm  for  construction  of  symbolic  sequence  and  recon¬ 
struction  of  chaotic  trajectory  from  symbolic  sequence  only  for  a  particular  set  of  parameters. 
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maps  for  weakly  chaoSc  mode  f%  222?Ts  , shows  these 
by  points  is  changing.  Even  the  coordinate  of  the 

are  permitted  ^  which^r^bhed^  analyzed  which  symbolic  sequences 

bolic  sequence  for  long  chaotic  tr^^ly  “^ly^ed  a  sym- 

reported  as  histograms  for  modes  with  F=0  233  and  V-O  99^  !?•  o-?  i. 
mstograms  demonstrating  which  seauence<?  nf  t,  '  ^  0.225.  Fig.  27  shows 

ics.  Sequences  are  numerated  by  natural  numbers  ®  ^  are  pemutted  by  generator  dynam- 
dex  of  a  sequence  is  simplylcfaS&m  correspondence  is  indicated  in  Table  2;  In- 

as  binary  nmber.  if  the  sequence  is  treated  as  a 


Table  2 


Number 


Sequence 
LL  (00) 
LR  (01) 
RL  go) 
RR  gi  ') 


As  one  c«t  see  in  Fig  27.  a  sequence  with  number  0,  i.e,  sequence  XT,  is  prohibited  in  these 

Ofpanicubr  ch^dc  oscillator  in  some  prac- 
elements  of  correspon^rMc  .Sc^^^r'ln  T'T’  ““  ‘"“.^tst  succeLn 

distorted  in  communication  ch  ’  ^  ^  ’  signal  or  a  ^gn^ 


0.6 


0,4 


0,2 
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(a) 


0.8  I — . r-  •■'.•■ 

oi— — 

00  01^ 


(b) 


Fig.  27.  Histograms  °f  ^'"quencfes  ^permitted  binary  symbolic  sequences  oflength 

2  m,(I4)  for  (a)  V  =  0.233  and  (b)  V  =  0.225.  ^ 


Fig.  28  Histograms  of  relative  frequencies  of  permitted  binary  symbolic  sequences  of  length 

4  in  (14)  for  (a)  V  =  0.233  and  (b)  V  =  0.225 . 


Ap^er  interesting  queston  concerns  the  possibility  of  discerning  different  modes  of  a  gen- 
erator  based  on  microtnode  by  means  of  analysis  of  symbolic  dynamics.  For  example,  do  the 
sets  ^permitted  sequences  differ  for  the  modes  mentioned  above?  To  answer  this  question  it 

of  Icngft  four.  Bu,  for  hotter  visualization  w^cSr  se' 
q  ences  of  length  five.  Histograms  for  these  sequences  are  presented  in  Fig.  28. 

uiLSors'^  wth  F  =  0.233  sequences  are  permitted  with 

numbers  30  md  31,  i.e.sequei.cesnM/i.  if^  H  iARii! ,  which  are  prohibited  for 

mode  TOth  F- 0.225.  Hus  difference  enables  to  discern  these  two  modes  in  the  oscillator 
simply-by  determining  prohibited  sequences  for  the  analyzed  signal. 


5.2.  Symbolic  description  for  TWT’System 

Construction  of  compact  information  description  for  TWT-system  (8)  is  a  much  more  com. 


Fust,  It  IS  necessary  to  reconstruct  .the  phase  space  of  the  system.  Second,  it  is  necessary  to 
find  a  region  of  modes  with_  relatively,  simple  dynamics.  Third,  even  for  low-dimensional 
modes  the  search  for  convement  Poincare  section,  which  could  give  good  unimodal  quasi- 
one-aimensional  map,  is  a  non-trivial  task. 


To  construct  the  succession  map  the  embedding  in  three-dimensional  space  is  used  This  is  an 
essental  restriction  of  our  method.  .  . 


The  first  method 

Here  an  example  of  construction  of  a  binary  symbolic  description  is  considered  for  several 
modes  with  fixed  parameters  a  =  2.05,  p  =  1.45.  Parameter  s  is  varied  in  the  range 
[0.4. ..0.45].  ® 

For  8=0.45  a  period  doubling  bifurcation  takes  place  (Fig.  29).  Fig.  30  shows  the  change  of 
quasi-one-dimensional  map  with  decreasing  of  s.  Starting  from  e  =0.447  the  map 

becomes  single-valued.  (Fig.  30,  31). 
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appearance  of  chaotic  oscillalions  Wtim  fh  a  -  ■  to  threshold  of 

succession  map^  Zltrdi !  “f “ '‘“^"“od,  the 

plex  parameterizStoToTSe  ShMuc&fof'^  f  ^  “ore  com- 

one-dimensional  map  to  eliminate  amhim.itv™  reversible  transformation  of  quasi- 

to  cornet  ate  genLting  partitiom  SuSn  Ts  if 

rithmofreconstructionofsa:S:;“is^^^^ 

conclusion  branch  that  favors  the 

Nevertheless  some  problems  described  bellow  can  be  considered  even  for  such  map. 

map'^Sfof  ^TO  mZtXiffiTrTchTaM  tb'^^  '-526  so  that  the 

Further  estimates  show  that  the  seniiA  ♦  ®  S*ueratmg  partition  can  be  introduced. 

the  initial  cha?S  “  “  a  close  connection  with 

time.  If  s^boIfseq^craSe  coScteff  th^*f  chaotic  trajectories  diverge  with 
coincide  but  from  a  trajectories,  then  first  several  symbols 

symboUc  sequences  coincide,  the  ctaareXXJitSe”' “““  “'“beta  of 


Fig.33.  Map  for  cross-section  in  Fig.  .32. 


Fig.  34-  Root-mean  square  error  of  initial  condition  identification  as  a  fimction  of  the  number 
of  coinciding  members  of  binary  sequences  belonging  to  different  trajectories  (a)  for  system 

(8)  and  (b)  for  system  (15). 


Therefore,  if  we  show  that  with  an  increase  of  the  number  of  coinciding  members  of  con¬ 
structed  bin^  sequences  the  chaotic  trajectories  become  closer  to  each  other,  we  would  get 
an  argument  (but  not  the  proof)  supporting  the  viewpoint  of  close  connection  between  the 
constructed  binary  sequence  and  initial  chaotic  trajectory. 

At  Fig.  34a,  horizontal  axis  contains  the  number  of  coinciding  members  of  binary  sequence, 
vertical  ms  -  discrepancy  in  initial  value  averaged  for  many  trajectories.  Here  the  fiajectories 
are  considered  with  the  specified  number  of  coinciding  members. 

As  one  can  see  from  Fig.  34a,  the  more  members  of  binary  sequences  coincide,  the  more  ac¬ 
curately  initial  points  of  the  trajectories  match. 
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'  P- im.  ™s  is 
formation  about  initial  conditions.  ^  ^  average, the  same  amount  of  in- 

stScS  5ls\tmbrS.^'"'  '»=  co„- 

“”=•  I-  -  of  iwr-syste.  ^.e  curve 

the  difference  betwS^MmiderS  nSoiTmci  mo  ^  Wtoximately  generating.  Nevertheless 
spond  to  the  chaotic  attractor  structure  could  ,  „  which  exactly  corre- 

reconstruction  of  a  signal  from  bina^-  s^bolic  sequ»«  is’mqm-rS'’' 

The  second  method 

The  essence  of  the  method  is  as  follows.  Denote  by  f,  the  moments  when  the  trajectorv  of  the 
system  moves  through  some  value  x-x  1  ^nott,  r^^ +•  •  .  , ,  ^^ajectory  of  the 

ingswedenoteasr  =r  -r  “‘'"'td  between  succession  cross- 

map  for  the  *®.  obtained  quasi-one-dimensional 

from  that  described  in  th  K  essentially  different 

®  ^^^^'^^^'d^ensional  maps  are  almost  the  same. 

theyrr”ne'Fta1hf  ^Zito  ^  P«fetable  compared  to 

be  more  stable  m  not.  ®  *'  “"^«ed  sequeZe  cZ 


Fig.  35.  Quasi-one-dimensional  map  r,.,(t,)  constructed  from  one-dimensional  realization 


sequences  of  lengft  fc^^e  „ed  in  dIS 

As  one  c^  see,  despite  that  methods  of  construction  of  binary  sequences  are  comnletelv  dif 

leneth  four  the  ^  sequences  among  sequences  of  length  two.  Among  sequencL  of 

length  four  the  following  sequences  are  prohibited;  LLLL ,  LLLR ,  RLLL . 


.  .  "  00  •  01  .  10  11 
Fig.  3  6.  Histogr^s  of  relative  frequencies  of  permitted  sequences  of  length  2  for  two  differ- 
ent  methods  of  construction  of  binary  sequence  of  the  system  (8). 


Fig.  37.  Histograms  of  relative  frequencies  of  permitted  sequences  of  length  4  for  two  differ¬ 
ent  methods  of  construction  of  binary  sequence  of  the  system  (8). 


6.  CONCLUSIONS 

The  conducted  investigation  consists  of  three  major  stag^^^ 

klystton.  MMC  and  other  atandard  RF  wSTngS®  fetf'  TWT, 

TO  and  low-dintensional  ohaosS  ean  g.” 

and“o”ttSfed“®Tof  has  shown  that  suggested  hitherto 

demonstrate  eompietcIcUhtSg  ^S 

and  most  fully  investigated  are  the  models  of  *  ■  e  most  covered  m  the  literature 

•  two-cavity  klystron  - 

•  generator  based  on  TWT 

•  microtriode  generator 

•  Pierce  diode  generator 

•  i^^^i^srd-wave  tube  generator 


vices’artd  can  ensmeTt^hXeSS^  design  of  the  de- 

for  qualitative  understanding  of  the  processes.  nsidered  genenc 

!!,!ns'^of°fk''“™t*‘r‘,?'®  knowledge  of  the  symbolic  description  of  the  trajectoiy  and  the  eoua 

«2t2:,r2=r  •'  *• 

investigations  was  conducted  to  construct  symbolic  description  for  the  model  of 


additional  reversible  map  t<  iminate  ambivalence  of  the  quasi-om  inensional  map  is  re- 
qmred.  Note  that  the  attractor  does  not  have  a  simple  structure  that  uniquely  identifies  the 
choice  of  the  succession  map.  Correlation  dimension  of  the  chaotic  attractor  for  which  the 
introduction  of  symbolic  dimension  is  possible  does  not  exceed  2.1.  There  is  no  success  in 
construction  of  symbolic  description  with  the  aid  of  quasi-one-dimensional  map  for  higher- 
dimensional  modes.  . 

The  results  obtained  state,  that  symbolic  description  in  the  investigated  systems  allows  (with 
some  limitahons).  identification  of  chaotic  modes  (for  example,  different  modes  of  the  same 
generator).  This  is  done  by  means  of  analysis  of  permitted  and  prohibited  sequences  of  a 
given  length.  Additionally,  uniqueness  of  the  set  of  permitted  and  prohibited  sequences  can  be 
used  in  various  applications  such  as  information  transmission  based  on  chaotic  synchroniza¬ 
tion  and  location  detection. 
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